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A new framework for the analysis of unstable oscillator flows is explored. In linear set¬ 
tings, temporally growing perturbations in a non-parallel flow represent unstable eigen- 
modes of the linear flow operator. In nonlinear settings, self-sustained periodic oscilla¬ 
tions of finite amplitude are commonly described as nonlinear global modes. In both 
cases the flow dynamics may be qualified as being endogenous, as opposed to the ex¬ 
ogenous behaviour of amplifier flows driven by external forcing. This paper introduces 
the endogeneity concept, a specific definition of the sensitivity of the global frequency 
and growth rate with respect to variations of the flow operator. The endogeneity, de¬ 
fined both in linear and nonlinear settings, characterizes the contribution of localized 
flow regions to the global eigendynamics. It is calculated in a simple manner as the local 
point-wise inner product between the time derivative of the direct flow state and an ad¬ 
joint mode. This study demonstrates for two canonical examples, the Ginzburg-Landau 
equation and the wake of a circular cylinder, how an analysis based on the endogeneity 
may be used for a physical discussion of the mechanisms that drive a global instability. 
The results are shown to be consistent with earlier ‘wavemaker’ definitions found in the 
literature, but the present formalism enables a more detailed discussion: a clear distinc¬ 
tion is made between oscillation frequency and growth rate, and individual contributions 
from the various terms of the flow operator can be isolated and separately discussed. In 
particular, in the context of nonlinear saturated oscillations in the cylinder wake, such an 
analysis allows to discriminate between the quasi-linear dynamics of fluctuations around 
a time-averaged mean flow on one hand and the effect of harmonic interactions on the 
other hand; the results elucidate why a linear analysis of the mean flow in this particular 
case provides accurate predictions of the nonlinear dynamics. 


1. Introduction 

Global instability in flows denotes the possibility of a spontaneous bifurcation from 
a steady flow state to a time-periodic state of synchronised oscillations in the entire 
flow field. A commonly observed scenario is that of a supercritical Hopf bifurcation, 
where linearly unstable perturbations of small amplitude first experience exponential 
growth, until nonlinear effects lead to amplitude saturation. The final time-periodic flow 
state is named a nonlinear global mode ( Huerre fc Monkewitz|[l9M ); the exponentially 
growing small-amplitude perturbations in the early stage of the bifurcation correspond to 
eigenmodes of the linearized flow operator, traditionally called linear global modes. The 
attribute global is used here to designate an analysis that resolves all non-homogeneous 
flow directions, as opposed to a local Ansatz, which implies the approximation of locally 
parallel flow. 

Nonlinear global modes are usually obtained as asymptotic oscillatory states from time- 
resolved numerical simulations, whereas linear global mode analysis requires the solution 
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of linear eigenvalue problems. Complex eigenvalues represent the temporal growth rate 
and the oscillation frequency; the associated eigenfunctions characterize the spatial distri¬ 
bution of fluctuation amplitude and phase. Linear global mode analysis is now routinely 
applied to two- and three-dimensional flow configurations. Beyond the primary question 
whether or not perturbations at small amplitude are unstable, a physical discussion of 
linear global modes is usually centered around two questions: (i) what are the physical 
mechanisms that give rise to unstable growth, and (ii) by what means can instability be 
reduced or enhanced? The first of these questions addresses the endogenous (or intrin¬ 
sic) flow behaviour, the second question concerns the control of those dynamics through 
exogenous (or extrinsic) manipulation. 

Huerre fc Monkewitz| (|1990 ) describe the conceptual notion of a ‘wavemaker’ (a word 


first used by Monkewitz [1990 ) as the region where instability waves are intrinsically 


generated in globally unstable flows. The interpretation by Koch (19851 of global in¬ 


stability in a wake already uses the same principal idea. Chomaz et al. (19911 derive a 


formal criterion for the global frequency selection in the context of the linear Ginzburg- 
Landau equation, based on the local absolute instability properties. Their formalism is 
rooted in a WKBJ approximation of instability wavepackets developing in a weakly non¬ 
parallel open flow. Within this approximation, local instability waves with upstream- and 
downstream-oriented group velocity emanate from a streamwise station, the ‘wavemaker’ 
location, where the two mode branches can be matched by means of a non-physical an¬ 
alytic continuation of the dispersion relation, defined as a function of a complex spatial 
x-coordinate. Such intrinsically generated waves grow and decay as they propagate. While 
the localized ‘wavemaker’ selects the frequency and drives the global instability mode, 
it is in general not characterized by large oscillation amplitudes. The spatial separation 
of the region where waves are generated and the region where they reach their maxi¬ 
mal amplitude is caused by convective instability mechanisms in a local sense, or by the 


non-normality of the linear Navier-Stokes operator in a global sense (Cossu & Chomaz 


1997 ). 

A quantitative theory of frequency selection in nonlinear systems, still based on the 
assumption of slow streamwise flow development and for the Ginzburg-Landau model 
equation, has been proposed by Couairon & Chomaz (1997) and by Pier et al. (1998). 


These studies draw on the theory of front dynamics (van Saarlos 1988, 1989), leading 


to the simple criterion that the nonlinear global mode frequency is given to first order 
by the absolute frequency at the upstream boundary of a (locally) absolutely unstable 
flow region of finite extent. Subsequent applications to wake flows ( |Pier fc Huerre|[2001 


|Pier|200^ Chomaz|2003 ) suggest that the accuracy of this criterion is only limited by the 
non-parallelism of the base flow over the distance of amplitude saturation. The transition 
from (upstream) convective to (downstream) absolute instability marks the ‘wavemaker’ 
location within the framework of this nonlinear model. 

Linear global modes in non-parallel flows may now be computed directly, without 
the need for the hypothesis of slow streamwise development. However, the notion of a 
cause-and-effect relation between different streamwise regions is lost along with this ap¬ 
proximation, and the localisation of a ‘wavemaker’ region within a global structure must 
be accomplished through new criteria. The sensitivity of the linear eigenvalue (frequency 
and growth rate) with respect to localized changes of the flow operator provides the 
appropriate concepts for a formal definition of a global ‘wavemaker’. Yet the sensitivity 
problem may be posed in several ways, depending on the physical premise of what ‘drives’ 
an instability. 


Giannetti & Luchini (2007) provide a discussion of the cylinder wake instability based 


on the structural sensitivity of the unstable linear eigenmode. The structural sensitivity. 
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in their definition, quantifies how an eigenvalue is affected by the introduction of localized 
forcing of a given perturbation quantity, proportional in strength to the same or another 
perturbation quantity. It thereby provides a measure in every point in space for the effect 


of internal feedback between perturbations. Giannetti & Luchini (20071 conjecture that 


those regions where an altered coupling induces the strongest change of the eigenvalue 
must also be the most significant regions for the action of internal feedback mechanisms 
that underpin the genuine eigenmode dynamics. At present, this formalism is arguably 
the most commonly accepted definition of the ‘wavemaker’ in a global analysis framework. 
The concept is quite naturally extended to nonlinear global modes by way of Floquet 
theory ( Luchini et a/.|[2008 |. However, one inconvenience of this approach is that it does 
not distinguish between frequency and growth rate, as the Cauchy-Schwarz theorem is 
invoked in order to define an upper bound for the drift of the modulus of the eigenvalue. 
Another stems from the large number of feedback relations between the various flow 
quantities that may prove to be significant. The formalism does not allow to single out 
the influence of specific terms in the flow equations. 


Marquet et al. { 2008) investigate the sensitivity of the linear cylinder wake instability 


with respect to localized modifications of the base flow. As far as linear instability is 
linked to the interaction between perturbations and the base flow, it may be argued that 
such an analysis is well suited to identify the principal flow regions where instability 
originates. The formalism distinguishes between frequency and growth rate, yet it is 
clearly cast in the form of a control problem. The question how an instability mode is 


affected by exogenous alterations, be it alterations of the internal feedback (Giannetti & 


Luchini 2007) or of the base flow (Marquet et al. 2008), is conceptually different from 


the question how its endogenous dynamics come into being. 

The objective of the present paper is to propose a variant of the sensitivity problem for 
linear as well as nonlinear global modes that identifies more directly those endogenous 
eigendynamics. For the sake of comparison and validation, these concepts are demon¬ 
strated for the two traditional test settings used in the literature on wavemakers: the 
one-dimensional Ginzburg-Landau equation and the two-dimensional cylinder wake. It 
is hoped that the formalism will be useful for the analysis of physical instability phenom¬ 
ena in a wide range of applications. 

The paper is organised as follows. Section 2 documents linear and nonlinear global 
mode results for the Ginzburg-Landau equation and for the cylinder wake. This section 
does not contain genuinely new results, but rather serves as a repertory and review. 
The configurations discussed here are used in the following as examples in order to 
demonstrate the proposed formalism. The endogeneity concept is introduced for linear 
settings in §3, and its application for the analysis of linear global modes is demonstrated 
for the two example configurations. The extension of the formalism to fully nonlinear 
situations is laid out in §4. Gonclusions are given in §5. An appendix addresses the 
implications of general inner products for the analysis. 


2. Linear, nonlinear, direct and adjoint global modes 

The evolution equation of a flow variable q(x, t) is considered in the general form 

Hdtq = A7(q), (2.1) 

where AC is a nonlinear operator, and B is an operator of very simple structure that 
indicates on which component of q the time derivative applies. In what follows, it will 
always be assumed that B is self-adjoint. 
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A base flow q;, is a steady solution of the nonlinear equation, JV{qi,) = 0. The linear 
stability of such a steady base flow is investigated by superposing small-amplitude time- 
dependent perturbations qb The dynamics of these perturbations is governed by the 
linear equation 

eatq' = /:q,q', (2.2) 

where the linear operator is obtained as the linearisation of Af around the base flow: 
Af{q_b + eq') = Af{(lb) + e£q^q'. Eigenvalues Xj and associated eigenfunctions of 

are obtained as solutions of the eigenvalue problem 


Xj B 


(2.3) 


The eigenmodes q) (x, t) = ^j(x) exp(Ajt) form a complete basis for the range of £qj. In 
physical terms, the real and imaginary parts of a complex eigenvalue represent temporal 
growth rate a and frequency w of a linear eigenmode. The convention A = ct — iw is 
adopted here, i.e. the frequency is given by the negative imaginary part of A. 

For the purpose of deriving an adjoint equation associated with the direct equation 


(2.2), the following inner product of vector-valued functions is introduced: 


{/(x,t),g(j 


.t)}= n n 

Jo Jn 


,t) • t/(x,t) dxdt, 


(2.4) 


where the star denotes the complex conjugate. The adjoint linear operator £1^^ is then 
derived by requiring 


{q^S(9tq' - £qtq'} = {-Bdtq< - £jj,^q'f,q'} . 
The adjoint equation associated with \2.2\ is thus found as 

-Bdtc^ 


(2.5) 

( 2 . 6 ) 


The linear adjoint operator has eigenmodes qj(x, t) = ^j(x) exp(A*t). Adjoint eigen¬ 
values A* are the complex conjugate values of the corresponding direct eigenvalues Xj. 
Adjoint eigenfunctions </>](x) satisfy 


A*S0]=£^^0]. (2.7) 

The direct and adjoint eigenfunctions form biorthogonal sets, 

(c/)],6c/)fe) (2.8) 

with respect to the purely spatial inner product 


(/(x),g(x))= / /*(x) •flr(x)dx. 


(2.9) 


In the definition of Huerre & Monkewitz (1990), a nonlinear global mode is a time- 


periodic solution of the nonlinear evolution equation (2.1), denoted here by qo(a;,t). In 
the examples considered in the present investigation, such a nonlinear global mode is the 
limit-cycle solution that is reached as the result of amplitude saturation of an initially 
growing linear eigenmode. The nonlinear global mode has a real-valued global frequency 
ujg = 27r/r, with T being the cycle period. 

The linear stability of this time-periodic solution is investigated by considering the 
temporal evolution of small perturbations q'(x, t), governed by the linear equation 


Bdtc{{x,t) = £q„(t) q'(a;,t). 


( 2 . 10 ) 
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The tangential linear operator Cqg(t) is obtained by linearizing the nonlinear operator Af 
around the time-periodic solution qo(t)- As Hqoit) is also time-periodic, a fundamental 


set of solutions to (2.10) is given by t) exp(^jt) (see looss & Joseph 1997). The 

Floquet modes (x, t) are T-periodic in time, and the associated Iloquet multipliers 
exp((*,T) characterize the temporal growth or decay of a mode over one cycle period. The 
real part of Q is in fact the Lyapunov exponent, the imaginary part of corresponds to 
a variation of the fundamental frequency. 

An interesting property is that the time derivative of the time-periodic solution, 9tqo, is 
a neutral Floquet mode of the autonomous periodic operator /Iqo(t). This property results 
from the phase invariance of the time-periodic solution. Both qo(x, t) and qo(x,t -|- 5t) 
represent T-periodic solutions of the nonlinear equation (2.1), therefore their difference 
= qo(x, t + 5t) — qo(x, t) is T-periodic as well. A Taylor expansion for small 5t gives 
Jq = i9tqo(x,It follows that 9tqo(x,t) is a T-periodic solution of (2.10), and is 
therefore a neutral Floquet mode: 


■0i(x,t) = atqo(x,f), Ci = 0- 

An adjoint tangential operator can also be defined by requiring 

{-Satq^= {q^ S^tq'- ^qp^q'} , 


from where follows the adjoint equation associated with (2.10), 

-eatq^(x,t) =/:|j^(j)q^(x,t). 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


The linear operator is also T-periodic. Equation (2.13) has fundamental solutions 

in the form of adjoint Floquet modes ■0j(x, t) exp(Ct). The sets_^ and xpf, are again 


biorthogonal, with respect to the spatio-temporal inner product (2.4), such that with a 
suitable normalisation they fulfill 




(2.14) 


Only the adjoint Floquet mode ip\ associated with the neutral Floquet mode and 
Cl = 0 will be used in the following analysis. 

2.1. Ginzhurg-Landau equation 

The Ginzburg-Landau equation has often served as a simple model for flow instability 


dynamics ( Huerre|2000 ). Its scalar state variable q{x, t) only depends on one spatial coor¬ 
dinate, yet its dispersion relation permits a double branch point that allows to distinguish 
between absolutely and convectively unstable situations, analogous to open shear flows. 
The nonlinear Ginzburg-Landau equation is used here in the form 


dtq = Af{q) = -Uda:q + n{x)q + jdxxq - I3\q\'^q ■ 


(2.15) 


The operator B, in the general notation (2.1), in this case is simply the identity, and the 


operator Af is given by the right-hand side in (2.15). It is composed of terms representing 


convection, linear reactive sources, diffusion, and cubic nonlinearity. Following |Gossu fc:| 
Chomaz (1997), constant convection and diffusion parameters are chosen, U = 6 and 


7 = 1 — i, while a parabolic variation of the reactive parameter is prescribed as 

yL{x) = fio + 0Afj,2X^ ■ 


(2.16) 


Different values of yiQ will be used in the following, while /i 2 = —0.1 is maintained 
throughout. This variation yields a strong local stability of the system far from x = 0. 
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All following calculations are performed on an interval x G [—40,40], and in all cases 
the fluctuation amplitudes are indeed negligibly small at the numerical boundaries. The 
nonlinearity parameter is chosen as /3 = 1 — i, except when linear situations /? = 0 are 
considered. 


2.1.1. Linear Ginzburg-Landau 

Note that the zero state g = 0 is a steady solution of the nonlinear Ginzburg-Landau 


equation (2.15). Linearisation around that state yields the linear Ginzburg-Landau equa¬ 
tion 

dtq = Lq = -Udxq + ix{x)q + ^d^xq , 


(2.17) 

which is equivalent to (2.15) with /3 = 0. A subscript 0 could be attached to £, denoting 


the zero base state, but this will be omitted in the following. The associated adjoint 
operator is obtained as 


-t 'y*dxxq^ ■ 


(2.18) 


As derived by Ghomaz et al. (1987), the leading eigenmodes of the direct problem 


X(j) = C(f> and of the adjoint problem are found as 


A — /To — 


47 V 

-H 2 I 

2 ’ 

(2.19) 

'Ux 

I-H 2 

( 2 . 20 ) 

.27 V 

27 2 J ’ 

' Ux 

/-M2 

( 2 . 21 ) 

. 27 * 

V 27 * 2 y ■ 


For fixed values of U, 7 and ^ 2 , the real-valued parameter completely determines the 
eigenvalue, but it has no influence on the eigenfunction shape. The direct and adjoint 
eigenfunctions are shown in hgurefl] Their m axima are located at a; = 8.6 and x = — 8 . 6 , 
respectively. According to Gossu & Ghomaz (1997), global instability arises when /rg > Hc, 
with a critical value 

/arg 7 ' 


Me — 




-b 


-M 27 

2 


cos 


V 2 


( 2 . 22 ) 


The local instability properties are given by [Ghomaz et al(| (1988). The model is locally 
stable wherever fi{x) < 0 , convectively unstable for 0 < fi{x) < and absolutely 

unstable for ^{x) > The extent of these regions depends on the value of ^g. 

2.1.2. Nonlinear Ginzburg-Landau 


Pier et al. (1998) report that the Hopf bifurcation in a nonlinear Ginzburg-Landau 


system is supercritical. Nonlinear global instability therefore follows from linear global 
instability, ^g > /ic, and the nonlinear global mode emerges after saturation of the 
growing linear global mode. A strongly supercritical setting is considered here, with 
Ho = ‘^IJ'c- 


The direct nonlinear global mode is computed by time-stepping (2.15), starting from 


the linear global mode shape ( 2 . 20 ) as an initial perturbation at low amplitude, until an 


asymptotic time-periodic state qo is reached. The linearisation of Equation (2.15) around 


the oscillating state qo is accomplished as detailed in Hwang (2015). The nonlinear term 
/3|qpg is properly linearized by augmenting the complex state variable with its complex 
conjugate, i.e. writing out the nonlinear and linear equations in terms of state vectors 
(g, g*)^. In the present study, the direct linear equation is in fact never needed, except 
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Figure 1. Direct and adjoint eigenfunctions (2.20 2.211 of the linear Ginzburg-Landau equa¬ 
tion. Absolute values are shown as thick blue lines, real parts as thin red lines. Solid lines 
represent the direct mode, dashed lines represent the adjoint mode. The scaling is such that 
= 1 . 


for the derivation of the adjoint linear equation, 

- dtq^ - Udxq^ - d^xq^ - 2/3*|go(Op9^ + P*ql{t)q^* = 0- 


(2.23) 


With the known time-periodic solution qoft), this equation is numerically integrated back¬ 
ward in time, until convergence towards a periodic solution gj is reached. The existence 
of such an asymptotic solution is guaranteed, because this is the adjoint Floquet mode 
associated with the neutral direct Floquet mode qi , defined as the time-derivative of the 
nonlinear time periodic solution qi = dtqo- Details on Floquet theory can be found in 


looss & Joseph (1997). 


A Crank-Nicolson scheme is used for the time integration of the nonlinear direct and 
the associated linear adjoint problem. The use of upwinding finite-difference stencils 
(seven points) for the spatial discretisation of the direct problem, and downwinding for 
the adjoint problem, is essential in order to achieve the required numerical accuracy. 

Figure]^ a) shows the nonlinear global mode go as it is recovered after transients have 
disappeared. The amplitude envelope has the emblematic shape of an ‘elephant’ mode 
( |Pier et aL||1998 ), with a sharp upstream wavefront and a softer downstream decay. 
According to WKBJ theory, the front should be situated at the upstream boundary of 
the absolutely unstable region in x. In the present case, absolute instability prevails in 
the interval x G [—10,10], and the foot of the front is indeed placed around x = —10. 
The adjoint Floquet mode, represented in figure [^&), has significant amplitudes only 
upstream of the direct wavepacket, with a maximum near x = —13. 

It is to be noted that the nonlinear global mode of the Ginzburg-Landau equation 
only contains one single frequency, because the nonlinear term is of such a form that 
it does not generate harmonics. As a result, the time signal in each point shows pure 
sinusoidal oscillations with zero mean. This property proves to be very convenient for 
all further discussion, whereas the harmonics and the non-zero mean oscillations that 
are characteristic for the nonlinear dynamics of the Navier-Stokes equations add further 
complexity to the analysis (see [|4.4[). 


2.2. Two-dimensional cylinder wake 

The wake of a cylinder is the most commonly used example of an oscillator-type shear 
flow. The critical Reynolds number for onset of self-sustained vortex shedding, Rcc = 47 


according to the experiments by Provansal et al. (1987), has been repeatedly recovered 
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Figure 2. a) Direct nonlinear global mode of the Ginzburg-Landau equation, with parameters 
as given in the text. Thick blue line: amplitude envelope; thin red line: snapshot of the real part 
of q. b) Associated neutral adjoint Floquet mode. Same line styles as in a. 


with high precision as the critical value for the onset of linear global instability (see 


for instance Barkley 2006 Sipp & Lebedev 2007). Furthermore, the nonlinear global 


instability of wakes is rather accurately predicted by local theory (Pier 2002 Chomaz 


2003). The cylinder wake has been chosen to illustrate the sensitivity studies by Giannetti 


& Luchini (2007), by Marquet et al. (2008), by Luchini et al. (2008) and by Luchini & 


Bottaro (2014). The case at Re = 80, well above the instability threshold, is considered 


in the present study. 

In the general notation (13, the state vector q = {u,p) now gathers the velocity 
vector and the pressure field, which together satisfy the incompressible Navier-Stokes 
equations, and the operators B and J\f are defined by 


B = 


AA(q) = 


-(u. 


V)u — Vp - 

Vu 


Re ^ Au 


(2.24) 


2.2.1. Linear instability of the cylinder wake 

The cylinder wake problem permits a non-trivial steady solution q;, = (uf,,pt,), sat¬ 
isfying Af{qb) = 0, which will serve as a base flow. Linear perturbations q' = (u',p') 


developing on this base flow are governed by the linear equations (2.2) with 

-(Ufc. V)u'-(u' 




V)u, 
V • u' 


Re Au'— Vp' 


(2.25) 


Eigenmodes of this linear operator are obtained numerically as described in |Sipp &:| 
Lebedev| ( |2007[ ). Only one unstable mode is found for Re = 80, with an eigenvalue 
X = a — loj = 0.1018 — 0.7852i . The streamwise velocity of the real part of this mode is 
shown in figure]^ a). The black line represents the stagnation-point streamline of the base 
flow, demarcating the recirculation region. The streamwise velocity of the associated 
adjoint eigenmode is displayed in figure [^&). As discussed in detail by Giannetti & 


Luchini (2007), Marquet et al. (2008) and others, the direct and adjoint eigenmode are 


localized upstream and downstream of the recirculation region, respectively. 


2.2.2. Nonlinear instability of the cylinder wake 

As the growing linear eigenmode reaches finite amplitude levels, the nonlinear terms 
become significant and the cylinder wake settles into a saturated periodically oscillat¬ 
ing state, the Benard-von Karman vortex street. This time-periodic solution qo(x,t) 
is obtained numerically by time-marching the nonlinear equations with a semi-implicit 
second-order temporal discretisation, and a spatial discretisation identical to the one used 
for the eigenvalue problem. At each temporal iteration, an unsteady Stokes problem is 
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(a) (b) 
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Figure 3. Unstable linear eigenmode of the cylinder wake at Re = 80. The real frequency is 
uJr = 0.785, and the temporal growth rate is Wi = 0.102. a) streamwise velocity perturbations 
u of the direct eigenmode; b) adjoint streamwise velocity perturbations of the associated 
adjoint mode.The black line indicates the recirculation region in the steady flow. 



Figure 4. Nonlinear global mode of the cylinder wake at Re = 80: a) total streamwise velocity 
U + w, h) adjoint streamwise velocity perturbations of the associated adjoint Floquet mode. 
Both snapshots are taken at the same instant. The black line indicates the recirculation region 
in the mean flow. 


solved with a preconditionned Uzawa algorithm (see Cuvelier et aZ.|1986 ) implemented in 
the FreeFem++ software. At Re = 80, the global frequency is found to be ojg = 0.9957, 
which is to be compared to the frequency of the linear eigenmode, w = 0.7852. The 
nonlinear correction to the global frequency is significant in this supercritical setting. A 
snapshot of the total streamwise velocity of the nonlinear global mode is shown in figure 
l^a), where the black line now represents the stagnation-point streamline contour of the 
time-averaged flow. In the same way as discussed in §2.1.2| for the nonlinear Ginzburg- 
Landau equation, see (2.13 2.23), an adjoint mode associated with the neutral Floquet 


mode of the nonlinear periodic state can be obtained by backward time-stepping of the 
adjoint tangential equation. Its streamwise velocity component vf is shown in figure|^6), 
taken at the same instant as the direct flow field in figure]^ a). 


3. Active flow regions in linear global modes 

3.1. Sensitivity of the eigenvalue 

The sensitivity of an eigenvalue measures how this value varies in response to changes 
of the operator. In the present context, only linear sensitivities are considered, i.e. all 


variations are assumed to be infinitesimally small. If the original linear equation (2.2) has 


eigenmodes that satisfy XjBcfj = Ccf>j, a small perturbation eSC of the operator leads to 
a perturbed eigenvalue problem, 


(Aj -l- eSXj)B(^(pj -|- e5(f)j) — (>C -I- e6£-){cj)^ -I- cScfA , 


(3.1) 


which at order e can be rearranged to give 


{XjB - C) S(f)g = -SXjBcfg + 5C cfg. (3.2) 

The Fredholm alternative states that this inhomogeneous problem in 5(f)j has a solution 
if and only if the right-hand side term is orthogonal to the nullspace of the left-hand side 
operator, i.e. (^j, —SXjBcfj +5C (f>j) = 0. If the direct and adjoint modes are normalized 
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such that ((/)J , B4>j) = 1, this condition leads to 

a, = (0j,<5£0,). 


(3.3) 


On the basis of (3.3), Giannetti & Luchini (2007) consider spatially localized perturba¬ 


tions of the operator, 


SC = (5(x- xo)Co , 


(3.4) 


where Cq represents some artificially added coupling between the various flow variables 
at the location Xq. If cf)j{xo) and (/)j(xo) are understood to be vectors containing the 
n flow variable values at Xq (n = 1 for the Ginzburg-Landau equation and n = 3 for 
the cylinder wake), then Cq is represented by an n x n matrix, Cq, and the eigenvalue 


variation is obtained from (3.3) as 


5 \, 


= cjC: 

Xo ^3 


Xxo) • Co • </>, (xo). 


(3.5) 


Taking the norm of Cq to be unity without loss of generality, application of the Cauchy- 
Schwarz theorem yields an upper bound for the modulus of the eigenvalue variation, 
induced by an operator variation at Xq ( [Giannetti fc Luchini|[20T7 ): 

(3.6) 


l<5Aj|xo ^ II^Uxo)|| ||</',(xo)|| 


Marquet et al. (2008) define the operator variation in (^3.3b specifically as being due to 


variations 5U of the base flow, SC — {VuC)SU. This definition allows to quantify how 
a given small modification of the base flow, localized or distributed, alters the frequency 
oj and the growth rate a. Both approaches represent the mathematical formulation of a 
well-posed question, based on different interpretations of what constitutes a ‘wavemaker’: 


in the case of 

Giannetti & Luchini 

2007 

), it is the localized ‘internal feedback’ between 

perturbations, whereas in the case o: 

Marquet et al. 

(2008 

) it is the feeding of perturbation 


growth on base flow energy. In both approaches, the answer is sought by probing the 
system with exogenous modihcations of the operator structure. 

3.2. Endogeneity analysis of linear global modes 
With the question in mind how a localized region in the flow contributes to the global 


dynamics, we note that the admittance of any arbitrary operator Cq in (3.4) may be 


too general for the purpose of identifying the specihc interactions that are inherent in 
the linear Navier-Stokes operator. Retaining the idea of considering the sensitivity of 
the eigenvalue with respect to localized changes of the operator, we stipulate that those 
changes preserve the local structure of the operator. This naturally leads to choosing 


SC = (5(x — Xq) C . 


(3.7) 


The variation of the operator at Xq is chosen to be proportional to the original operator 
itself in that same location. The sensitivity with respect to such variations quantifies 
directly how much the eigendynamics in a given point in space contribute to the frequency 
and to the growth rate; it is therefore suitable for an investigation of the endogenous 
global dynamics. We call this specific sensitivity 


E{x) = </>]*(x) • {Ccf>j) (x) 


(3.8) 


the endogeneity of the eigenmode (Aj, 4>j). Its computation is straightforward if the direct 


and adjoint eigenmodes as well as the operator are available. The dot-product in (3.8) 


again only denotes the scalar multiplication of two vectors containing the various state 
variables in one point in space. 




























Identifying the active flow regions that drive linear and nonlinear instabilities 11 




X X 

Figure 5. Endogeneity distribution of the leading linear Ginzburg-Landau eigenmode, (a) 
real part related to growth rate, (fe) negative imaginary part related to frequency. 


An essential property of the endogeneity is that its integral over x is equal to the 
eigenvalue Aj-: 

[ A(x)dx = [ cfY ■ dx = [cf] , Ccf^) = {cf \, = A,. (3.9) 

Jn Jq 

It is important to note that, while any spatial distribution can be normalised to yield 
any integral scalar value, the endogeneity is the unique quantity that represents local 
contributions to the eigenvalue. For instance, the endogeneity allows to exclude any 
point in space from the integration, and the result reflects how the eigenvalue is altered 
due to the missing local contribution. 

An endogeneity-based analysis clearly distinguishes between the promotion of unstable 
growth, contained in the real part of £'(x), and the frequency selection, given by the 
negative imaginary part. For the sake of clarity, let 


Act(x) = 3? [i?(x)] and A,^(x) = —9 [Fl(x)] 


(3.10) 


be defined, such that A(x) = — ii?^(x), analogous to A = cr — iw. The distinction 

between these two components is of great importance for a physical discussion, and the 
following examples will show that E„ and Ei^ in general present quite different spatial 
structures. Furthermore, with the definition (3.8) it is straightforward to decompose the 
operator £, for instance into convection, diffusion and other terms, and to examine the 
individual contributions of these separate parts. Such a decomposition will be discussed 
in §3.4|for the cylinder wake. 


3.3. Example 1: linear Ginzburg-Landau equation 
The endogeneity formalism is first applied to the global instability modes shown in figure 
with the parameters as given in § 2.1.1 Choosing = Be as defined by (|2.22), the 


system is marginally unstable in a global sense. The endogeneity is found by multiply¬ 
ing (ff* with C(j) in every point x, where the eigenfunctions 4){x) and (/>^(x) are given 
analytically by (2.20 2.21), and the operator L is written out in ( 2.17[ ). 

Both parts of the endogeneity, E^ and E^, are shown in hgure|^ both are even func¬ 
tions, with their maximum values at x = 0. The integral of E(x) is exactly equal to 
the associated eigenvalue, A = —4.398i. In the E^{x) distribution, shown in figurea), 
negative and positive regions exactly counterbalance each other, totalling a zero growth 
rate. The largest contribution to both the frequency selection and the growth rate stems 
from the region around x = 0, not from regions where the magnitude of either the direct 
or adjoint eigenmodes is large (compare to figure [^. 

This result is compared to the saddle point criterion given by Chomaz et al. (1991). In 
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Figure 6. Endogeneity distribution of the linear cylinder wake instability, a) i?o-(x), related to 
growth rate; 6) Eui{x.), related to frequency. The spatial integral of Ecr equals the growth rate 
<j — 0.102, the spatial integral of E,^ equals the frequency uj = 0.785. 


their model, the ‘wavemaker’ location is defined by a saddle point of the local absolute 
frequency ujq in the complex x plane, d^uJo = 0. In the present case, one finds ujo(x) = 
i^{x) — zC/^/dy, and therefore dxioo = idxfi, with a single saddle point precisely at x = 0. 


The ‘wavemaker’ location according to the criterion of Chomaz et al. (19911 is identical 


with the maximally endogenous location in this example. But whereas their WKBJ- 
based criterion identifies a singular location, the endogeneity quantifies contributions to 
the growth rate and the frequency from any point in the domain, thereby characterizing 
a distributed ‘wavemaker’. 

3.4. Example 2: unstable linear global mode of the cylinder wake 
The endogeneity of the unstable linear eigenmode 4> of the cylinder wake, displayed in 
figure [^a), is now examined. It is readily computed by point-wise multiplication of the 
complex conjugate of the adjoint eigenmode 4>', displayed in figure [^b), with Ccf) where 
the definition of C is given in (2.251. As the divergence of perturbation velocity is zero 


in every point, the continuity equation along with the adjoint pressure vanishes from 
the endogeneity definition, and one is left with 

E{x) = — u'^* • [(ub • V)u] — u'^* • [(u • V)u;,] — u^* • Vp -b Re~ 


.-lut* 


Au, 


(3.11) 


where again it is understood that the left-hand side is to be evaluated in every point 
X, such that all vectors only contain two scalar elements (x- and y-components). The 
endogeneity of the linear cylinder-wake instability mode is displayed in figure The 
distribution of ifcr(x), depicted in figure [^a), shows where the temporal growth rate is 
generated, whereas i?^(x), shown in figure [^b), indicates how the various flow regions 
influence the global frequency selection. 

A few general conclusions can be inferred from figure Mainly, it is observed that 
the endogeneity (real and imaginary parts) is concentrated around the shear layers of 
the separation region, delimited in the figures by black lines. The frequency selection is 
clearly concentrated in two symmetric maximum locations. This part of the endogeneity 


resembles the quantity displayed by Giannetti & Luchini (2007) in their figure 17. The 
distribution of Ea-{'x.) however bears a more faceted structure. Some regions are positive, 
contributing to global instability, others are negative, thus stabilising the eigenmode. The 
entire flow downstream of the separation region has practically no influence on frequency 
and growth rate, consistent with the conclusions of Giannetti fc Luchinij (2007) and 
Marquet et al. (2008). 


A more insightful analysis of physical mechanisms may be based on a decomposition 


of the endogeneity into contributions from the various left-hand side terms in (3.11). In 


the given order, these terms account for the effects of base flow convection; production 
through base flow shear; pressure forces and diffusion. Their spatial distributions (real 
parts only, reflecting contributions to the growth rate) are shown in figure]^ The effect 
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(a) (b) 




Figure 7. Real parts of the individual endogeneity terms in ( 3.11|| (from left to right), a) 
Convection by the base flow; b) production through base flow shear; c) pressure forces; d) 


diffusion. The sum of these contributions gives Eafil), shown in £gure|6 
is used here. 


The same color scale 


of convection by the base flow is dominantly stabilising (figure]^). In the outer vicinity 
of the separation bubble, the downstream convection of perturbations counteracts their 
capacity of in situ growth, and renders the instability more convective. Inside the recir¬ 
culation region, the upstream convection has the opposite effect. The production term 
(figure [^) provides the principal source of global unstable growth. 

The workings of pressure forces are not as obvious to interpret in physical terms, 
and their integrated net contribution to and is exactly zero in an incompressible 
setting. Since the adjoint velocity field is divergence-free, V • ul = 0, it is easily found 
that 

f • Vp dx = — f (V-u^)pdx = 0 (3-12) 

Jfi Jq 

Yet this term contributes strongly to the local values of the overall endogeneity. Tenta¬ 
tively, it may be argued that the role of the perturbation pressure gradient is to enforce 
the continuity condition, thereby causing a perturbation volume flux across the shear 
layer. As the instability perturbations tend to shorten the length of the separation bubble 
(manifest in the nonlinear mean flow), the stagnation-point streamline is forced toward 
the symmetry line, lessening the convective effect outside the bubble, and enhancing it 
inside. 

The diffusion term is globally stabilizing, although inside the separation bubble near 
the cylinder it provokes a weak destabilisation. This seems to be the consequence of 
viscous transport of perturbation velocity into the shear layer. However low the ampli¬ 
tude of the diffusive contribution, it is crucial for accurately determining the instability 
threshold. The Reynolds number is expected to have two distinct effects on the insta¬ 
bility. First, the stabilizing effect of the perturbation diffusion should weaken when the 
Reynolds number is increased. This is confirmed in figure|^a), which shows the trends of 
the convection, production and diffusion contributions to the growth rate, as functions 
of the Reynolds number. The stabilizing diffusion effect lessens with increasing Reynolds 
number. Secondly, the steady base flow is influenced by the Reynolds number, affecting 
the instability mechanism via the production and convection terms. Individually, these 
contributions appear to be dominant, but with opposite effects on the growth rate. As 
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Figure 8. Net contributions of individual terms in the endogeneity definition (|3.11|), (a) to the 
linear growth rate, and (6) to the global frequency, as a function of the Reynolds number Re. 
Convection (-C-), production (-p-) and diffusion (-d-) contributions are shown, as well as the 
sum of convection and production (c + p) and the total values of a and uj (thick red lines). The 
pressure gradient contribution is zero for all Reynolds numbers. 


the Reynolds number is increased, the destabilisation by the production term overcomes 
the stabilisation by the convection term. The combined effect of these competing terms, 
represented by a line marked c + p in figure |^a) , is then comparable in strength to the 
diffusion effect. On the other hand, the role of the diffusion in the frequency selection 
is all but negligible, as seen in figure [^b). As for the growth rate, the dominant con¬ 
tributors are the production and convection terms. Both contribute here to increase the 
frequency. Their combined effect (c -|- p) gives a good prediction of the linear frequency 
(red line), especially at high Reynolds number. 


4. Active flow regions in nonlinear global modes 

4.1. Sensitivity of the global frequency 


The general sensitivity formulation for nonlinear time-periodic oscillations of a bifurcated 
flow state adopted here is similar to the analysis by Luchini et al. (20081, with some vari¬ 
ations in the notation. Consider a nonlinear global mode qo(x,t), time-periodic solution 
of (2.1), with fundamental frequency uig. In order to make the frequency explicitly visible 
in the equation, the time variable is rescaled as r = Wgt, such that qo(x, r) is 27r-periodic 
in T and satisfies 


Wg d^qo = A/'(qo) . (4.1) 

Just as in the linear case of §3.1[ small variations of the left-hand side operator cause 
variations of the solution, including its frequency: 


(Wg + e6ojg)B dr {<10 + eJqo) = (AA -b eJA/')(qo + eJqo) . (4.2) 


With the introduction of the tangential linear operator /Iqo(t) and its neutral Floquet 
mode ipi = drqo, defined in ( 2.10[[24T| ), variations are governed at order e by the relation 

{ujgB - /Iqo(t)) S^pl = - 6 ujg B1P1+ 6JV{qo) , (4.3) 


The Fredholm alternative can be written out with the aid of the neutral adjoint Floquet 
mode V’l (see normalised to give {^p\,B^f}l} = 1, which leads to 

6ujg = {t/>|,(5A/'(qo)} • 


(4.4) 
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Figure 9. Endogeneity E^i (thick blue line) of the nonlinear global mode of the Ginzburg— 
Landau equation, as shown in figure!^ com pared wi t h the quantity h{x) (thin red line), which 
is obtained with the formalism described by Hwang (20151 according to (4.7l. The dotted line 
marks the location Xca, and the dashed line traces the envelope shape of the nonlinear global 
mode. 


If 'ip\ is known, then the impact of any small operator variation SJV on the global fre¬ 


quency can be immediately evaluated from (4.4). In practice, t/ij is obtained in the fol¬ 


lowing way: first, the nonlinear equation is numerically integrated by time-stepping until 
the periodic nonlinear global mode regime is fully attained. Then the (linear) adjoint tan¬ 
gential equation is stepped backwards in time, starting from an arbitrary initialization, 
over as many cycles of the nonlinear global mode as necessary. During this backward-in¬ 
time integration, the adjoint solution converges asymptotically towards the sought-after 
neutral mode. 


4.2. Endogeneity analysis of nonlinear global modes 
Analogously to the linear case in § |3.2[ the endogeneity of a nonlinear global mode is 
defined by considering variations of the nonlinear operator that preserve its structure but 
are localized in space, and also in time, since the nonlinear operator is time-dependent, 


SAf = (5(x — xo)(5(t — To)Af. 


(4.5) 


The inner product in (4.4) involves integration in x as well as in t. Integration over the 


Dirac functions yields the expression for the influence of spatio-temporal variations in 
the operator on the frequency selection. 


E^{x,t) = ^ ■0|*(x,t) • A/'(qo(x,r)) 


(4.6) 


Only the frequency selection in nonlinear global modes is considered at present, because 
it is assumed here that both quantities, in practice, are obtained from a flow solver in 
the form of real-valued variables. The scaled time r denotes the temporal phase within 
an oscillation cycle. 


4.3. Example 1: the nonlinear Ginzburg-Landau equation 
The nonlinear global mode and its neutral adjoint Floquet mode of the Ginzburg-Landau 
equation (figure have been discussed in (2.1.2 The associated endogeneity, according 
to ( |4.6| ) , is shown in figure It is noted that the endogeneity is independent of r in the 
present case of the Ginzburg-Landau equation. 

The Ei^ distribution is nonzero in a narrow region around x = —9, and its integral in 
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X and over one period 0 ^ t ^ 27r is exactly equal to the global frequency, uig = 3.70. 
According to Pier et al. (19981, the global frequency of a nonlinear Ginzburg-Landau 
system with slowly varying coefficients is selected at the location Xca, where the local 
instability changes from upstream convective to downstream absolute. In the limit of 
marginal global instability, the global frequency is then predicted to correspond to the 
absolute frequency at Xca- In the present example, one finds Xca = —10, and the absolute 
frequency at this location is ujQ(xca) = 4.50. The chosen parameter configuration (see 
^2.1) is strongly globally unstable, and the parameter /i varies significantly in x around 
therefore the present case does not respect the limiting assumptions of [Pier et al.\ 


(1998), and the frequency prediction is quite inaccurate as a result. The endogeneity 


however provides a clear and accurate picture of the frequency selection process. The 
maximum contribution to the global frequency is found at x = —9, indeed not far from 
the ‘wavemaker’ location Xca as defined by Pier et al. (1998). 


The present results may be compared to the structural sensitivity of nonlinear global 
modes as defined by Hwang (2015), who adapted the formulation of Luchini et al. (|2008 ) 
to analyze nonlinear global modes of the Ginzburg-Landau equation. According to Hwang| 
(2015), worst-case variations of the global frequency due to added ‘closed-loop pertur¬ 
bations’ (synonymous to ‘internal feedback’), written in the notation of this paper, are 
characterized by the spatial distribution 


h{x) = 


q^{x) tplix) 


N 


T /*oo 


with N = 


{-ipl* qo - tpl qo) dx dt. 


(4.7) 


0 ^ —oo 


Figure ^ compares this distribution, multiplied with 2ujg for consistent scaling, with 
the endogeneity The two curves are very similar, and in particular the position of 


their maxima is identical. Again it is pointed out that the analysis of Hwang (2015), in 


contrast to the endogeneity formalism, considers modifications of the operator that do 
not preserve its original structure. 


In conclusion, just as in the linear analysis of (3.2 the results obtained from the non¬ 
linear endogeneity analysis are consistent with the ‘wavemaker’ definition from classical 
asymptotic theory, but they go further in a quantitative description of the dynamics, be¬ 
cause the endogeneity fully accounts for the effects of non-parallelism and supercriticality. 
It is also consistent with recent formulations that describe the structural sensitivity, but 
it reveals the endogenous dynamics that are specific to the operator under consideration. 

4.4. Example 2: Nonlinear frequency selection in the cylinder wake 
The endogeneity of the time-periodic flows developing in the wake of the circular cylinder, 
described in § 2 . 2 . 2 [ is explicitly obtained as 

Euj{y^, r) = uj • [-(uo • V)uo - Vpo + i?e"^ Auq] , (4.8) 

where qo(x, r) = (uo,po) is the 27r-periodic nonlinear solution and u| is the velocity 
component of the 27r-periodic solution of the linear adjoint equation. 

The spatio-temporal endogeneity (4.8) is integrated over one oscillation cycle, and 
the result, shown in figure [T^ demonstrates that the endogenous region resembles a 
front, localized around x = 2. Similarly to the earlier observations in the context of the 
Ginzburg-Landau equation, the endogenous region of the nonlinear global mode in the 
cylinder wake is located further upstream than that of its linear counterpart (compare 
to figure |^) . The separation line of the recirculation regions both in the base flow and 
in the time-averaged mean flow are depicted by grey and black lines, respectively, in 
figure [T^ Interestingly, the endogenous region of the nonlinear global mode appears to 
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Figure 10. Time-integrated endogeneity Eui{x.) of the nonlinear global mode in the cylinder 
wake at Re = 80. The nonlinear frequency is uig = 0.995. The black (resp. grey) lines indicate 
the recirculation region in the time-averaged mean flow (resp. base flow). 


be supported by the recirculation region in the mean flow, as opposed to the base flow. 
This observation suggests that further investigation of the nonlinear dynamics should 
focus on the analysis of fluctuations around the mean flow state. 


To this aim, the nonlinear global mode state is first decomposed as qo = q + q' 
into a time-averaged component q = {u,p)"^ and a zero-mean fluctuation com ponen t 
q' = (u',p')^. Introducing this decomposition into the Navier-Stokes operator (2.24), 
and time-averaging over one oscillation cycle, one obtains the steady nonlinear equations 
that govern the mean flow. 


{u ■ V)m -I- Vp — Re ^Au = —{u' ■ V)m', 
V • M = 0 . 


(4.9) 

(4.10) 


By subtraction, the unsteady nonlinear equations for the fluctuations are obtained as 


, dru' + (u ■ 'V)u' + (u' ■ V)u + \7p' — Re ^ Au' = —{u' ■ V)m^ -I- {u' ■ ^)u', (4.11) 


V • m' = 0. 


(4.12) 


The right-hand side forcing terms both in (4.9) and in (4.11) arise from the nonlinear 


interaction of fluctuations. The temporal mean of («'• V)m' forces the mean flow, whereas 
its zero-mean fluctuation forces the flow fluctuations. At small amplitudes of u', these 
terms are negligible, such that the steady and unsteady equations reduce to the base flow 
and linear perturbation equations. By contrast, when a linear eigenmode experiences 
exponential growth and reaches finite amplitude levels, the effect of the forcing terms 


in (4.9) and (4.11) will become significant. The right-hand side in (4.9) drives the base 


flow towards the mean flow. The right-hand side in (4.11) modifies the dynamics of 
the fundamental oscillations at frequency Wg, and it generates harmonic components at 
integer multiples of uig. Thus the effect of nonlinearity on the frequency selection can be 
attributed to two distinct origins. The first is a nonlinear deformation of the mean flow. 


which in turn modifies the left-hand side linear operator in (4.11). The second one is the 


nonlinear interaction of fluctuation harmonics, induced by the right-hand side forcing 


term in (4.111. 


The endogeneity definition (4.81 is now expanded in terms of mean flow and fluctuation 
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Figure 11. Contributions to the time-integrated endoge neity of the nonlinear global mode, (a) 
Quasi-linear dynamics around the m ean fl ow defined by (4.151 and (b) int erac tion between high 
harmonical components defined by (4.151. The color bar shown in figure 10 is also used here. 


The black line indicates the recirculating flow region in the mean flow. Re = 80. 

components. With this decomposition, as developed above, one obtains 
t) = u'^ • — {u ■ V)u — Vp + Re~^Au — (u' • V)u' 

— (u • V)u' — (u' • V)u — Vp' + Re~^Au' 


(4.13) 


- (u' • V)u' -b (u' • V)u' 

The first line in the above equation vanishes because the expression in square brack¬ 


ets is the momentum equation for the mean flow (4.9 By further splitting the adjoint 


velocity into mean and fluctuation components, denoted and respectively, the 
time-integrated endogeneity simplifies to 

^27r p‘2'T^ 

/ dr = / Em{x,T)dT+ / Eh{x,T) dr , 

Jo Jo Jo 


(4.14) 


with 


Ejn{x, t) = u^' ■ (— (u • V)u' — (u' • V)u — Vp' -b Re 
Eh{x,T) = u'^' ■ (-{W ■ V)u' -b (u' • V)u'^ , 


^Au') 


(4.15) 

(4.16) 


All contributions involving the mean adjoint velocity vanish. The component E^{x,t) 
highlights spatial regions where the quasi-linear dynamics around the mean flow con¬ 
tribute to the frequency selection, whereas Eh{x,T) identifies regions where the interac¬ 
tion of harmonic components influences the global frequency. 

The time-integrated components Em and Eh, for the nonlinear global mode in the 
cylinder wake at Re = 80, are shown in figure [m The sum of these two contributions, 
according to (4.14), gives the time-integrated frequency endogeneity displayed in figure 


El In the present configuration, the frequency selection process is clearly dominated 
by the quasi-linear dynamics around the mean flow, whereas harmonic interactions con¬ 
tributes only very weakly. The total contributions 

p p2.'K p p‘2'k 

Wm = / / Em{y^,T) dr dx, U}h= Eh{x,T)dTdx, (4.17) 

Jn Jo Jn Jo 

are reported in table for various values of the Reynolds number. For all values of the 
Reynolds number investigated in this study, the contribution uJm of the quasi-linear mean 
flow dynamics to the global frequency is greater than 97%. 
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Re 


Wg 

Wm(%) 

0Jh{%) 

50 

0.7676 

0.8159 

0.8087(99.1) 

0.0072(0.9) 

75 

0.7742 

0.9726 

0.9543(98.1) 

0.0182(1.9) 

80 

0.7852 

0.9957 

0.9768(98.1) 

0.0189(1.9) 

100 

0.7553 

1.0703 

1.0493(98.0) 

0.0210(2.0) 

125 

0.7246 

1.1382 

1.1127(97.7) 

0.0255(2.3) 

150 

0.6896 

1.1922 

1.1596(97.2) 

0.0325(2.8) 


Table 1 . Endogeneity-based decomposition of the nonlinear global mode frequency, for various 
values of the Reynolds number Re. Linear eigenmode frequency oj of the base flow; nonlinear 
global mode frequency ujg\ contributions to ujg from qu asi-lin ear dynamics in the mean flow 
(tUm) and from harmonic interaction {u)h), as defined by (4.171. 


This result explains a posteriori why a global stability analysis of the time-averaged 
flow accurately predicts the frequency of the time-periodic vortex shedding in the cylin¬ 
der wake (Barkley 2006). Mantic-Lugo et al. (2014) recently proposed a self-consistent 
nonlinear model based on the marginal stability of the mean flow, which accurately deter¬ 
mines both the mean flow and the frequency of the vortex shedding in the cylinder wake. 
However, the influence of higher harmonics in the frequency selection of nonlinear insta¬ 
bility is not negligible in all circumstances. Indeed, Turton et al. (2015) recently showed 
that in thermosolutal convection driven by opposite thermal and solutal gradients, oscil¬ 
lation frequencies of travelling convection waves can be predicted from stability analysis 
of the mean flow, but not those of standing waves. The application of endogeneity anal¬ 
ysis to the thermosolutal convection problem is not attempted here, but it is expected 
to reveal a stronger influence of the harmonic interactions in the frequency selection of 
standing waves. 

Finally, a decomposition of the endogeneity component (4.15) similar to the demon¬ 
stration in §3.4| characterizes the contributions of the various terms to the quasi-linear 


dynamics. For Re = 80, the contributions of production and mean flow convection to the 
global frequency are displayed in figures [T^ , b, and the contributions of diffusion and 
pressure gradient are shown in figures [T^ , d. The production contribution is dominant, 
concentrated in the shear layers of the mean recirculation region. Inside this region, the 
action of the pressure gradient further increases the frequency. Table summarizes all 
total contributions, of the production, convection and diffusion terms at 

various Reynolds numbers . The net integral of the pressure gradient term is always iden¬ 
tically zero. At all Reynolds numbers, the effects of diffusion and mean flow convection 
approximately compensate each other, and their balance is small compared to the strong 
contribution of the production term. 


5. Conclusions 

A novel sensitivity formalism has been introduced, named the endogeneity, which al¬ 
lows to precisely quantify the influence of each point in the flow field on the global fre¬ 
quency selection and on the promotion of unstable growth, both in the context of linear 
temporal eigenmodes and of nonlinear global modes. Its application has been demon¬ 
strated for the Ginzburg-Landau equation and for the wake of a circular cylinder, in 
linear as well as nonlinear settings. The results obtained have been shown to be con¬ 
sistent with earlier ‘wavemaker’ definitions, in particular the WKBJ-based saddle point 


criterion of Chomaz et al. (1991) and the structural sensitivities defined by Giannetti h 
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Figure 12. Spatial distribution of various contributions to the endogeneity of the quasi-linear 
mean flow dynamics: (a) production by the mean flow, (b) convection by the mean flow, (c) 
diffusion and (d) pressure gradient. The black line delimits the recirculation region in the mean 
flow. 


7?e I I ^tr. 


U) 


d 

m 


50 

0.8087 

0.9997 

0.0023 -0.1933 

75 

0.9543 

0.9435 

0.1567 -0.1459 

80 

0.9768 

0.9482 

0.1679 -0.1394 

100 

1.0493 

0.9999 

0.1696 -0.1203 

125 

1.1127 

1.0989 

0.1224 -0.1087 

150 

1.1596 

1.1721 

0.0853 -0.0978 

175 

1.1935 

1.2220 

0.0604 -0.0889 


Table 2. Decomposition of the quasi-linear contribution to the fresquency into three contribu¬ 
tions: the production term convection term u>^ and diffusion term in the quasi-linear 
mean flow dynamics. 


Luchini (20071, Luchini et al. (20081 and Hwang (20151. The novel aspect with respect 


to the latter sensitivity approaches arises from the specific form of operator variations 
that are considered: the endogeneity characterizes the sensitivity of the eigenvalue with 
respect to localized operator variations that preserve the specific structure of the original 
operator. This sensitivity may therefore be interpreted as the local contribution of any 
point in the flow field to the global eigendynamics, in terms of frequency selection and 
unstable growth. Contributions to these two parts of the eigenvalue are clearly distin¬ 
guished, contained separately in the real and imaginary parts of the endogeneity. Further 
analysis of the role of individual terms of the operator follows naturally within this frame¬ 
work. In particular, a decomposition of a nonlinear flow operator into time-averaged and 
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fluctuating parts gives insight into the role of quasi-linear dynamics developing on the 
mean flow versus the nonlinear interaction of harmonic fluctuation components. This 
latter part of the analysis, exemplified for the case of the cylinder wake, has important 
implications for the characterisation of nonlinear time-periodic flow states, in relation to 
recent investigations (Mantic-Lugo et ^|2014 Turton et qZ.||201'5). 
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Appendix A. Influence of the inner product 


The question may arise whether and how the choice of an inner product different from 
(2.9) affects the endogeneity definition. In particular, many flow problems are investi¬ 
gated in cylindrical coordinates, where the standard inner product includes the radial 
coordinate r as part of the volume element. For practical purposes, the procedure is 
outlined here for discrete eigenvalue problems, with a generalized inner product. 

Let 4>j be a discrete representation of the linear eigenfunction (pj, in a discrete space 


where the inner product {pi, 02) between two flow states is expressed as (pi Q(i> 2 - The 
matrix Q is typically diagonal, with real elements that represent volume elements of the 
mesh, and possibly any further weight functions. It will only be assumed here that Q is 
invertible, but even this condition can be relaxed with some additional effort. 

The discrete direct and adjoint eigenvalue problems are defined by 


—lUJjBcpj = L<pj , 


(Al) 
(A 2) 


The matrix Q on both sides of (A2) can be omitted. It is found that (pj = Q^cp^j is 
the adjoint eigenvector satisfying 




f H jA 


(A3) 


Let ^j{xk) be defined as the vector = J(x — Xk)L<pj for convenient writing, 

where x^ denotes the discrete mesh points. The endogeneity definition, written in terms 
of discrete vectors and operators, is then 


-iA(xfc) = {p],6{x- Xk) Cpj) = pl^Q^jixk) = 0)"$j(xfe). 




(A4) 

It follows from ( |A4[ ) that the endogeneity is invariant with respect to the choice of 
the inner product. This holds true also for the example of cylindrical coordinates: it 
seems unnecessary to account for the r factor. Without the need to specify a Q matrix, 
E[xk) may be computed directly from the discrete adjoint eigenvector, obtained from 
the transpose conjugate problem (A3). This will generally be the simplest option. 
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